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ABSTRACT 

We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early 
O supergiant f Pup with models that include the effects of porosity in the stellar wind. We explore the effects of 
porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented 
parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic 
clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify 
the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line 
shapes, and we show that only modest increases in the mass-loss rate (<40%) are allowed if moderate porosity 
effects (/loo < 1C) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled 
out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to 
both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate 
calibration for bright, nearby O stars. 
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1. INTRODUCTION 

Radiatively driven winds in massive stars have been known to 
be structured and inhomogeneous virtually since their discovery 
(Morton 1967), both from an observational and theoretical 
standpoint (Lucy & Solomon 1970; Brucato 1971). The same 
instability in the radiative driving force that is responsible for 
the inhomogeneous nature of these winds is also responsible 
for generating shocks resulting in X-ray emission (Lucy & 
White 1980; Owocki et al. 1988; Feldmeier et al. 1997b). Since 
these phenomena are so intimately related, it is no surprise that 
the question of the effects of wind inhomogeneities on X-ray 
radiative transfer has arisen as well. 

A good understanding of X-ray line profile formation in 
massive star winds is of great importance, as the degree of 
asymmetry in X-ray line profiles can provide an independent 
diagnostic of stellar mass loss (Cohen et al. 2010). From 
an observational perspective, the traditional mass-loss rate 
diagnostics in the optical, radio, and UV bands are easier to use, 
but all these diagnostics are subject to modeling uncertainties 
resulting from wind inhomogeneities. Ha recombination lines 
and radio free-free emission are both the result of collisional 
processes that scale with the square of density, and thus are 
sensitive to the degree of clumping in the wind. Ha emission 
originates near the stellar photosphere, while radio free-free 
emission originates farther out in the wind, and each of these 
measurements depends on the local degree of inhomogeneity 
(e.g., Puls et al. 2006). UV P Cygni absorption features probe the 
ion column density in the wind, but the derivation of a mass-loss 
rate from an ion column density requires accurate knowledge 
of the ionization balance (Massa et al. 2003; Fullerton et al. 
2006). Furthermore, inhomogeneities in the distribution of wind 
material in line-of-sight velocity complicate the construction 


of model absorption profiles (Owocki 2008; Oskinova et al. 
2007; Sundqvist et al. 2010, 2011). If X-ray line profiles are 
demonstrated to be relatively free of systematic effects from 
wind inhomogeneities, then X-ray spectroscopy may be used to 
calibrate the effects of inhomogeneities on traditional mass-loss 
rate diagnostics. 

Because X-ray radiative transfer in stellar winds is generally 
only affected by continuum photoelectric absorption in the cool 
bulk of the wind, the only possible way for inhomogeneities 
to have an effect is if the clumps are optically thick in the 
continuum and have a large interclump mean free path, leading 
to porosity. On the other hand, optically thin clumps will only 
affect diagnostics with opacities that scale with the square of the 
density, and so do not affect X-ray absorption. The goal of this 
paper is to empirically constrain the potential effect of porosity 
on X-ray line diagnostics. 

Feldmeier et al. (2003) and Oskinova et al. (2004) have 
developed a formalism for X-ray radiative transfer in a porous 
stellar wind, focusing on the case of flat clumps oriented parallel 
to the stellar surface. They find that porous winds can produce 
significantly more symmetric line profiles for a given mass-loss 
rate. Since the assumed mass-loss rate also influences the degree 
of profile asymmetry, the mass-loss rate and degree of porosity 
might be expected to become degenerate parameters of line 
profile fitting. This would reduce the utility of X-ray line shapes 
as mass-loss rate diagnostics, unless some other observational 
constraints could be placed on porosity lengths. 

In Owocki & Cohen (2006) and Sundqvist et al. (2012; here- 
after Paper I), we have developed a generalized formalism for 
X-ray radiative transfer in a clumpy stellar wind. In partic- 
ular, Paper I gives a unified treatment of spherical and flat 
clumps, and explores the effects of a distribution of clump scales. 
Paper I shows that porosity in models with either fiat or 
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spherical clumps tends to produce more symmetric line shapes 
for a given mass-loss rate. However, flat clumps result in profiles 
with a distinct shape, with a bump arising near line center (the 
“Venetian blind effect;” e.g., columns 3 and 4 of Figure 3 in 
Paper I). 

Oskinova et al. (2006) have applied their X-ray radiative 
transfer formalism to the calculation of emission line profiles, 
and have compared their model to the observed line profiles 
of four O stars observed with Chandra. They found that a 
model including porosity from flat clumps could qualitatively 
reproduce the observed profile shapes. However, they have 
not formally evaluated the quality of their fits, and have not 
quantitatively explored the tradeoffs between porosity and mass 
loss rate. 

In this paper, we aim to confront models of X-ray emission 
line profiles with high quality observational data to evaluate 
how well porous and non-porous models can explain the X-ray 
emission of O stars, to evaluate the effects of both flat and 
spherical clumps, and to quantify the tradeoffs between porosity 
and mass-loss rate in modeling the observations. We focus 
on the high quality archival observations of the bright, well- 
studied O star £ Pup, which has been observed by Chandra, 
as well as by XMM-Newton for nearly one million seconds. 
In addition to the high quality of the available data, £ Pup 
has the advantage of being a single, non-magnetic O star, and 
it has visibly asymmetric X-ray line profiles, indicating that 
absorption of X-rays by the wind is significant and measurable. 

This paper is organized as follows. Section 2 briefly reca- 
pitulates our X-ray emission line profile models. In Section 3 
we describe the observations, and our data reduction and fitting 
procedures. Section 4 gives the results of our line profile mod- 
eling. In Section 5 we compare our results to previous work. In 
Section 6 we discuss our results and give our conclusions. 

2. MODEL 

The models used in this paper are described in Owocki & 
Cohen (2001), with elaborations for porosity in Owocki & 
Cohen (2006) and Paper I. Here we give a brief recapitulation 
of these models and a description of their implementation. 

The emergent X-ray spectrum from the wind of a massive star 
is calculated assuming the X-rays are formed in a small fraction 
of the wind, while the cool bulk of the wind absorbs the X-rays 
as they propagate through it. The emergent luminosity at a given 
wavelength X is 

Lx = 4?r J dVr]xe ~ T , (1) 

where the integral is taken over the X-ray emitting volume of 
the wind, rjx gives the X-ray emissivity, and r is the X-ray 
optical depth due to bound-free transitions in the cool bulk of 
the wind, evaluated by an integral along the line of sight. The 
X-ray emitting volume is typically taken to have a lower radial 
cutoff Rq. In most of this work, we assume that there is no 
upper radial cutoff R m ax , but in Section 5 we instead assume a 
finite value for R mdx in order to compare our work with that of 
Oskinova et al. (2006). 

It is possible to incorporate porosity within the framework of 
Owocki & Cohen (2001) by making appropriate modifications 
to the optical depth term in Equation (1). In this paper, we 
use the implementations of isotropic and anisotropic porosity 
described in Paper I. The bound-free absorption of an X-ray 


photon emitted at z e in direction z is given by the optical depth 

poo 

r = Xeff dz, (2) 

where the effective opacity Xeff accounts for any porosity, /cir is 
calculated using the “inverse” (or Rosseland) law for bridging 
the optically thin and optically thick regimes 


X 1 + Tel ’ 

where x is the atomic opacity per unit length, proportional to 
the mass-loss rate and here characterized by the fiducial optical 
depth r* = Mk /(Ait R t v x ), with a mass absorption coefficient 
k and a wind terminal speed Voo- The clump optical depth is 

Td = xVImI> (4) 

where /i is the directional cosine for a photon impacting a 
clump; isotropic porosity is recovered by replacing /x with unity 
(Paper I). Note from Equation (3) that the effective opacity 
approaches the atomic opacity when r c i <gC 1 . For the porosity 
length h we assume the “velocity stretch” form hir^/h^ = 
v(r)/v oo, where the velocity field is given by the standard 
“/1-law” v(r)/Voo = (1 — R*/r)P . Thus, to evaluate X-ray line 
profiles accounting for porosity, it is only necessary to specify 
one additional parameter beyond r* and Rq : the terminal porosity 
length hoc. 

It is also possible to account for the effects of U V photoexcita- 
tion of metastable levels in He-like ions and resonance scattering 
of strong X-ray lines within this framework (Leutenegger et al. 
2006, 2007). Because the focus of this paper is to derive ob- 
servational constraints on porosity and its effects on mass-loss 
rate determinations, we do not use fits to He-like triplets for this 
purpose since they are blended. However, in our exploration 
of the effects of a finite upper radial cutoff R max in Section 5, 
we fit He-like triplet complexes to obtain constraints on R max . 
Also, with the exception of our fits to He-like complexes, we 
also do not account for resonance scattering, which might have 
a moderate effect on derived mass-loss rates, but will not affect 
our conclusions about the suitability of porous models or the 
quantitative tradeoff between porosity and mass-loss rate. 

All of the models are implemented as local models in the 
X-ray spectral fitting package XSPEC (Arnaud 1996). The code 
is freely available for download 6 and is issued under the General 
Public License. 

3. OBSERVATIONS, DATA REDUCTION, 

AND MODEL FITTING 

£ Pup was observed by Chandra on 2000 March 28-29 using 
the High Energy Transmission Grating Spectrometer (HETGS; 
ObsID 640, 67.7 ks effective exposure time). £ Pup was also 
observed several times in the last ~10 years by XMM-Newton, 
with a total effective exposure time in excess of 0.8 Ms; the 
list of observations used is given in Table 1. The Medium 
Energy Grating (MEG), High Energy Grating (HEG), and 
XMM-Newton Reflection Grating Spectrometer (RGS) spectra 
are presented in Figure 1 . 

The HETGS data were reprocessed using CIAO 4.2 and 
CALDB 4.3.1 in order to apply optimized contemporary charge 


6 http://heasarc.nasa.gov/xanadu/xspec/models/windprof.html 
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Table 1 

List of XMM-Newton Observations of f Pup with Net Exposure Times 


ObsID 

Date 

Exposure Time (ks) 
RGSl RGS2 

0095810301 

2000 Jun 8 

52.6 

51.0 

0095810401 

2000 Oct 15 

39.9 

38.5 

0157160401 

2002 Nov 10 

41.6 

40.2 

0157160501 

2002 Nov 17 

38.7 

38.7 

0157160901 

2002 Nov 24 

43.5 

43.5 

0157161101 

2002 Dec 15 

27.8 

27.8 

0159360101 

2003 May 30 

66.3 

66.3 

0163360101 

2003 Dec 6 

41.5 

41.5 

0159360301 

2004 Apr 12 

27.7 

27.7 

0159360401 

2004 Nov 14 

57.3 

57.3 

0159360501 

2005 Apr 16 

34.6 

34.6 

0159360701 

2005 Oct 15 

23.5 

23.4 

0159360901 

2005 Dec 3 

48.3 

48.2 

0159361101 

2006 Apr 17 

42.5 

42.4 

0414400101 

2007 Apr 9 

58.5 

58.5 

0159361301 

2008 Oct 13 

61.2 

61.2 

0561380101 

2009 Nov 3 

64.1 

64.2 

0561380201 

2010 Oct 7 

76.7 

76.7 

Total 


846.2 

841.8 


transfer inefficiency (CTI) calibrations. Standard procedures 
were used to produce spectra, response matrices, and ancillary 
response files. 

The RGS (den Herder et al. 2001) data were processed using 
the XMM SAS 11.0.0. rgsproc was used with the following 
significant modifications from default settings: pixel-by-pixel 
corrections to CTI based on contemporary diagnostic exposures 
were applied; periods of high background were rejected; events 
on CCD node interfaces, near the edges of CCDs, and next to bad 
pixels were not rejected; “cool” pixels that have been flagged 
as having consistently high CTI were rejected. The resulting 
spectra were then coadded using rgscombine. 

RGS spectra are affected by small, systematic shifts in the 
wavelength scale. The shifts are different for each observation, 
but constant within each observation. The shifts have an rms 
deviation of er = 8 mA (den Herder et al. 2001). The most 
recent calibration implements a correction to the RGS boresight 
which zeroes out the mean systematic wavelength scale shift 
in the data used for calibration (Coia & Pollock 2008). Since 
we have combined a large number of observations of £ Pup, 
the overall systematic wavelength scale shift is reduced, and the 
combined data set is broadened. For N independent observations 
of equal length, the systematic shift in the combined data set 
is expected to be of the order of 8 mA /y/~N. Since we are 
using 18 observations of roughly equal length, we estimate a 
net systematic uncertainty in the wavelength scale of ~2 mA. 

Based on the fitting results we report in Section 4, we 
believe that the corrections to the RGS boresight implemented 
in the most recent version of the SAS still result in relative 
differences between the RGS 1 and RGS2 wavelength scales of 
~3-5 mA. Thus, any quantitative measurement of line profile 
asymmetry using RGS data will still be dominated by the 
systematic wavelength scale uncertainty, and joint fits to data 
from both RGS instruments will have poor fit statistics which 
are dominated by the systematic disagreement in the relative 
calibration of these instruments. 

Regardless of these systematic uncertainties in the absolute 
and relative wavelength scales of the RGS instruments, the data 


are of extremely high statistical quality and are well suited 
to answering our two key questions: first, do the data show 
a bump near line center indicating that anisotropic porosity is 
important? Second, what is the quantitative tradeoff between 
porosity effects and mass-loss rate reductions? 

We fit the strong emission lines in the HETGS and RGS 
spectra with the line profile models described in Section 2 using 
XSPEC 12.6 (Arnaud 1996). Each line was fit separately rather 
than performing a global fit. The continuum was modeled locally 
for each line using a power law in energy with index 2, which is 
a good approximation to the shape of the continuum for a small 
range in energy. For the Chandra data the continuum strength 
was determined for a given spectral region by direct fitting while 
excluding all lines. For the RGS data this is not possible because 
of the extended wings of the instrumental line- spread function 
(LSF), so the continuum strength was fit simultaneously with the 
line profiles. For lines where data from both RGS and Chandra 
are usable, we fit them separately and compare results. For lines 
where data from both RGS instruments is usable, they are fit 
jointly, using only negative first order data. In all fits to HETGS 
data, we used only first order data, fitting positive and negative 
orders simultaneously. We fit MEG and HEG simultaneously, 
when HEG data were available. When possible, we fit the line 
over a wavelength range about 20% bigger than the Doppler shift 
associated with the wind terminal velocity of v^c = 2250 km s 1 
(Haser 1995). When there is a nearby line which could be 
blended with the line in question, we eliminate part of the wing 
of the line from the fit, a procedure which has been shown to 
have little effect on the fit and the derived parameters (Cohen 
et al. 2010). We explicitly note these cases in our discussion of 
individual lines at the end of this section. 

Closely spaced doublets in the Lya lines are fit with a single 
profile model centered at the emissivity-weighted wavelength 
of the two components. 

We assess the goodness of our model fitting and assign 
confidence limits using the C statistic (Cash 1979) for the low- 
background, low-count-per-bin Chandra spectra, and use the 
X 2 statistic with Churazov weighting (Churazov et al. 1996) for 
the higher background and higher count RGS data. Confidence 
limits are assigned using the A/ 2 formalism described in 
Chapter 15.6 of Numerical Recipes (Press et al. 2007). We 
report formal 68.3% confidence limits throughout this paper 
unless otherwise noted. 

Because of the systematic wavelength offset between RGS1 
and RGS2, the fits to lines where data were available for both 
RGS instruments usually had X„ > 2. This also occurred in 
some cases where models provided very poor fits to the data. 
XSPEC deems a model with X„ 2 >2 a poor fit and does not 
provide confidence intervals for model parameters, so we do 
not give them in our results for the two lines where even the 
best-fit models do not provide formally good fits. For models 
with large porosity lengths, many other lines also do not have 
formal confidence limits quoted. 

Although it is possible to extract meaningful constraints on 
model parameters from fits to blended lines, we have chosen 
to focus on a few strong lines with minimal blending in order 
to emphasize constraints on porosity and quantitative tradeoffs 
between porosity and mass-loss rate. Thus, we do not include 
fits to He-like triplets in our study of the effects of porosity, 
nor do we include the N vii Ly« line, which is blended with 
Nvi He/3. 

Below, we give notes on several of the lines which we fit. 
Also, note that for the RGS there are several lines which show 
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variations in the effective area due to hot pixels (visible as narrow 
“dips” in both the spectrum and model). These variations are 
correctly accounted for by the standard RGS calibration data 
and pipeline processing. We also tested the effect of excluding 
these bins from our fits, and found that our results did not change 
significantly. 

Mg xii Lya. The RGS has relatively poor resolving power 
at this wavelength and does not provide strong constraints on 
profile shape, so we only fit the HETGS data. 

Nex Lya. The Nex Lya doublet is blended with a nearby 
Fe XVII line at 12.124 A, which has a comparable temperature- 
dependent emissivity. It is so close in wavelength to the Ne 
X lines, which are themselves very close together, that the 
feature can safely be modeled as a single profile. In contrast, 
there are two additional iron lines — another Fe XVII line at 
12.266 A and an Fe XXI line at 12.284 A, which forms at higher 
temperatures — that may be blended with the red wing of the Ne 
X line. We therefore exclude the longest wavelength portion of 
this line from our fitting for both the Chandra and XMM-Newton 
data. 

Fe xvn 15.014 A. The RGS data for this line are moderately 
affected by blending from the wings of neighboring lines, 
especially the strong Fexvn line at 15.261 A, so we have 
excluded the red wing of this line for the RGS data sets. There are 
also a number of weak lines from the Rydberg series of O vm 
in this spectral region, as well as satellite lines from Fexvi. 
Neither of these are strong enough to have a significant effect 
on the profile shape. 

Fexvn 16.780 A. This line is slightly contaminated by the 
Rydberg series of Ovu, but this does not have a significant 
effect on the line profile shape. We use only the MEG data and 
not the RGS data, because of contamination from the extended 
wings of the RGS LSF from the nearby lines at 17.051 and 
17.096 A. 

O viii Lya. There is slight blending with the Rydberg series 
of N vi, but these lines are sufficiently weak that they do not 
have a significant effect on the O vm profile shape. Most of the 
MEG counts come from the negative first order, which falls on 
one of the back illuminated ACIS CCD chips, which has much 
higher detection efficiency. 

N vn Ly 1 6 . This line is somewhat weak, but valuable because it 
is unblended. The Chandra data are too weak to give meaningful 
constraints. Most of the RGS observations put part of this line 
on a gap between CCD chips of RGS 1 , but enough observations 
exist with small pointing offsets to give acceptable exposure 
over the whole line. This line falls on a failed CCD of RGS2, so 
no data are available from that spectrometer. 

4. RESULTS 

We fit each line in the XMM-Newton and Chandra grating 
spectra according to the methodology described in Section 3. 
As discussed in that section, we exclude from our analysis 
those lines that are too weak or too blended to provide strong 
constraints on the parameters of the line profile model, described 
in Section 2. For those lines that do provide useful results, 
Tables 2 and 3 summarize the results of the line profile fits that 
do not include any porosity. Tables 4-7 summarize the results 
of line profile fits including porosity. 

In this section, we give single parameter confidence limits for 
all parameters of interest. These single parameter confidence 
limits are determined with all other parameters of interest free. 
They are evaluated for a change in the fit statistic of 1.0 
(68.3% confidence, corresponding to la). These confidence 


limits are based only on the change in fit statistic under 
variation of model parameters and do not reflect systematic 
uncertainties. 

4.1. Comparison of RGS and Chandra Results 

We first fit each line without including the effects of porosity. 
As can be seen by inspecting the r* and Rq parameters and 
their confidence limits in Tables 2 and 3, we find very good 
agreement for the three lines common to both RGS and Chandra 
datasets. As an illustrative case. Figure 2 shows the fits to 
the Fexvn line at 15.014 A in the RGS and HETGS spectra. 
The fitted t* values are consistent with a mass-loss rate of 
~3 x 1 0” 6 Mq yr -1 (or a correspondingly lower mass-loss rate 
if higher wind opacity is adopted), as derived from the ensemble 
of Chandra HETGS lines in Cohen et al. (2010), and the fitted 
Rq values of ~ 1 .5 R» are consistent with the predictions of line- 
driven instability (LDI) simulations of wind shock formation 
(Feldmeier et al. 1997b; Runacres & Owocki 2002). The overall 
agreement between the r* values found in lines from these two 
datasets has recently been confirmed by Naze et al. (2012). 

4.2. Constraints on Porosity Models from Individual Lines 

The fits presented in the previous subsection show that the 
simple line-profile model without porosity provides good fits 
to both the Chandra and XMM-Newton data. The derived 
parameters are physically reasonable in the context of embedded 
wind shocks. However, it is possible that models that include 
the effects of porosity could also fit the data. We now turn 
our attention to line profile models that include anisotropic or 
isotropic porosity. 

First, we report on the fitting of models that assume 
anisotropic porosity, from radially oriented, flattened clumps. 
We first performed fits that allowed the terminal porosity length 
to be a free parameter. These fits generally prefer h 0 0 — 0 (i.e., 
no porosity, recovering the same results as in Section 4.1). To 
further explore the effects of anisotropic porosity on the line 
profiles, we next fit models with fixed values of the terminal 
porosity length, h a 0 — 0.5, 1, 2, and 5 R*. We find increas- 
ingly poor fits as the value of the terminal porosity length in- 
creases. To illustrate these trends, Figure 3 shows the Fexvn 
line at 15.014 A, with the best-fit models assuming each of the 
four values of hoc- Plots of fits to the other lines in the spec- 
trum of f Pup using anisotropic porosity models are shown in 
the Appendix. 

Tables 4 and 5 summarize the results of fitting anisotropic 
porosity models to all of the lines. To facilitate comparison with 
the non-porous model fitting results, the table repeats those 
results in the rows with h x = 0. For models that provide 
very poor fits (i.e., rejected with high probability compared 
to the best-fitting h a 0 — 0 models), we do not list formal 
uncertainties on the free model parameters, as the formalism 
for assigning these uncertainties assumes that the best-fitting 
models do indeed provide formally good fits. 

One clear result that can readily be seen in these tables is that 
in no case is an anisotropic porous model statistically preferred 
over the non-porous model. Even for those lines with a lower 
signal-to-noise ratio, the mildest anisotropic porosity models 
(with hoo = 0.5 R*) can be rejected at the 68.3% confidence 
level. 

Let us next consider models with isotropic porosity. In this 
case, there is some degeneracy between porous profiles with 
lower optical depths and non-porous profiles with higher optical 
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Figure 1. The extracted coadded positive and negative first order MEG (top panel) and HEG (second panel) spectra along with the negative first order RGS1 (third 
panel) and RGS2 (bottom panel) spectra (coadded 18 pointings). 


Table 2 

Wind Profile Model without Porosity: Fit Results to the Chandra HETGS Spectra 


Ion 

X 

(A) 

T* 

Ro 

(7?*) 

Normalization 
(10 -4 photons cm -2 s -1 ) 

C 

M)ins 

Mg xn Lya 
Ne x Lya 
Fe xvii 
Fe xvii 
O viii Fy a 

8.4210 

12.1339 

15.014 

16.780 

18.969 

1.22+-™ 

2- 01 - 2 24 
1.94+' 3 3 2 3 

3- OHfo 
3.00+- 34 

1.34-S 

1.55+- 12 

1 01 + ' 59 
I.U1_ oi 

1.22+' 37 

0.294 4 ™ 

2- 7F°o 9 9 

5 ' 24 - 2 l 4 7 

2.45+- 1 1 3 7 

3 - 70 -.35 

186.5 

191.4 

280.8 

174.9 

150.9 

188 

176 

308 

308 

130 


depths, as shown in Figure 3 in Paper I. We find that models 
with isotropic porosity provide better fits than the corresponding 
anisotropic porosity profile models. However, when we allow 
/too to be a free parameter, non-porous (h a 0 = 0) models are 
almost always preferred, and in no case is li 0 0 = 0 rejected at 
greater than 90% confidence for any one line. 


We then fit models with fixed values of terminal porosity 
length, hoc = 0.5, 1, 2, and 5 R *. Figure 4 shows fits of the 
isotropic porosity models to the Fexvn line at 15.014 A, and 
Tables 6 and 7 present fits to all the lines. Plots of fits to the other 
lines in the spectrum of £ Pup using isotropic porosity models 
are shown in Figures 12, 14, 16, 18, and 20 in the Appendix. 
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Figure 2. The Chandra MEG (left panel) and XMM-Newton RGS1 (right panel) data for the Fe xvn line at 15.014 A, along with the respective best-fit models. The 
vertical dashed line in each figure represents the laboratory rest wavelength of the line. The vertical dotted lines denote the Doppler shift associated with the wind 
terminal velocity of Voa = 2250 kins 2 Although the left panel shows only the MEG data, the HEG data were fit simultaneously. The right panel shows only the 
RGS 1 data; however, both RGS datasets were fit simultaneously. 

(A color version of this figure is available in the online journal.) 


Table 3 

Wind Profile Model without Porosity: Fit Results to the XMM-Newton RGS Spectra 


Ion 

X 

T* 

Ro 

Normalization 

X 2 

^bins 


(A) 


(R«) 

(10 -4 photons cm -2 s -1 ) 



Ne x Lya 

12.1339 

1.8L^ 

1 61 +15 
1 — .18 

3 in+ 08 
- 3 - 1 - .08 

24.0 

19 

Fe xvii 

15.014 

1.77 

1.57 

6.39 

129.5 

52 

O viii Lya 

18.969 

3.14 

1.01 

4.66 

204.3 

72 

N vii Ly fi 

20.910 

4.932“ 

1 41+- 62 
1 — .40 

1-662% 

40.5 

36 


Figure 5 shows confidence limit contours for the fits to the 
15.014 A line in T t -h oc parameter space, illustrating the tradeoff 
between porosity and wind optical depth. Modest porosity 
(/too < R*) is not strongly ruled out, but it does not have a strong 
effect on r*, and thus on mass-loss rate estimates from X-ray 
line profiles. 

Figure 6 shows the best fit models for the Fexvn 15.014 A 
line in the Chandra data without accounting for the instrument 
response. This allows evaluation of the level of degeneracy 
of hoc and r* in producing roughly comparable line profiles. 
Even for a small porosity length of hoc = 0.5 Ah, the anisotropic 
porosity model shows a significant change in shape near line 
center, with the bump becoming more pronounced for larger 
porosity lengths. On the other hand, the isotropic porosity 
models with small h a 0 do not look very different from the non- 
porous model. However, the shape does become increasingly 
different as porosity length increases. The changes in profile 
shape support the same conclusions as Figure 5: for isotropic 
porosity, small porosity lengths can produce similar line shapes 
that adequately fit the data, while large porosity lengths are ruled 
out; and for anisotropic porosity, even small porosity lengths 
change the profile shape so that it does not fit the data. 

4.3. Effects of Porosity on Inferred Mass-loss Rates 

In Section 4.2 we have shown for several lines in the spectrum 
of £ Pup that anisotropic porosity models do not fit the data 
for any significant porosity length, and that isotropic porosity 
models do not fit for large porosity lengths. However, isotropic 
porosity models with moderate porosity lengths are formally 
allowed for some individual lines. 


Cohen et al. (2010) have fit the ensemble of measured r* 
values for a non-porous model of £ Pup by using the mass-loss 
rate as an adjustable derived parameter, together with a model 
wind opacity. We can apply the same technique to the r* values 
measured for different porosity lengths to obtain an estimate of 
X-ray derived mass-loss rate as a function of assumed porosity 
length. Together with constraints on the porosity length and 
geometry, whether from X-ray profile modeling or from other 
observations, this relation gives a quantitative constraint on 
the contribution of porosity to uncertainties in mass-loss rates 
derived from fitting X-ray line profiles. 

To this end, it is most instructive to consider not the mass-loss 
rate itself, but the fractional change in mass-loss rate relative to 
a non-porous model. The derived mass-loss rate scales linearly 
with measured r*, and thus the fractional change in derived 
mass-loss rate must equal the fractional change in measured r*. 
Thus, as a simple way to estimate the fractional change in mass- 
loss rate, we calculate the fractional change in measured r* for 
the ensemble of lines for a given assumed porosity length. The 
estimated mass-loss rate increase for isotropic porosity models 
is roughly linear in hoc, with A M/M ~ 0.4 for h 0 0 = R *. 
Large porosity lengths indeed produce large changes in inferred 
mass-loss rate compared to non-porous models, and anisotropic 
porosity models have a stronger effect than isotropic models. 
However, all of the models with strong effects are ruled out by 
our fit results. 

5. COMPARISON WITH PREVIOUS RESULTS 

One of the main results of Section 4 is that anisotropic 
porosity models are disfavored, as well as isotropic porosity 
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Wavelength (A) Wavelength (A) 

Figure 3. The Chandra MEG (left column) and XMM-Newton RGS (right column) measurements of the Fexvn line at 15.014 A, identical to the data shown in 
Figure 2, but with best-fit anisotropic porosity models superimposed. The characteristic “Venetian blind” bump at line center (see Figure 4 of Oskinova et al. 2004 
or Figure 3 of Paper I) can be seen in all these models, with the feature getting stronger with increasing h^. The effect of porosity on the profile shapes can be seen 
qualitatively even for the models with the smallest porosity lengths, and even with the moderate resolution of the RGS data. These anisotropic porosity profile models 
do not provide a good fit to the data. 

(A color version of this figure is available in the online journal.) 
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Table 4 

Wind Profile Model with Anisotropic Porosity: Fit Results to the Chandra HETGS Spectra 


Ion 

X 

^OO 

T* 

Ro 

Normalization 

C 

Abins 


(A) 

(«*) 


(R*) 

( 10~ 4 photons cm ~ 2 s"" 1 ) 



Mg xn Lya? 

8.4210 

0 

1 22 +-78 
i.ZZ _44 

1-34*1 

0.294+° 0 2 2 4 2 

186.5 

188 



0.5 

1 72 + - 77 

A - 7 Z -. 8 i 

1 27 + - 17 
1 7 —.14 

0-281+° 0 5 0 ° 8 

189.8 

188 



1 

1 46+1-58 
1 .60 

1 45+-14 
1 — .16 

0 2 OS + 018 

U.jUj_ 033 

191.5 

188 



2 

1 ,83+ 2 95 8 

L49+ 13 , 

0.298+°o 23 4 

194.1 

188 



5 

9 97+6.19 
z - 3/ -1.48 

1-54*1 

0.297+^3 

198.6 

188 

Nex Lya? 

12.1339 

0 

2 01 + ‘ 27 
Z - U1 — .24 

l-45 4l 0 3 8 

2 71 + 09 
z - 7 a -.09 

191.4 

176 



0.5 

2.92+^j 

1 49+09 
a -^-.09 

7 79+.09 
Z ' /z -.09 

208.8 

176 



1 

3-60+i 4 ! 

1 ^/^+-08 
l-5o_ os 

2 71 +10 
Z ‘ 7 1 —.08 

223.8 

176 



2 

500-1.01 

1 65 +07 
1-0 -.09 

2 72+ 10 
z - 7 —.09 

245.0 

176 



5 

11.32+” 4 

i 77 + 0.11 
17 7 - 0.10 

2 76 +010 
z - 7 —0.09 

277.6 

176 

Fe xvii 

15.014 

0 

L94 +- 32 

1-55*'. 12 

5.24+1 

280.8 

308 



0.5 

2 S 1 + - 68 
z -' ?1 -.40 


5-23+1 

293.9 

308 



1 

3.02 +' 84 

1 65 +10 
A -° -.10 

5-22 + 23 

303.8 

308 



2 

• 3 - /o - 1.00 

1.74+ 0 ! 9 ! 

5 22 +-23 
* 5-ZZ — .18 

317.4 

308 



5 

6-95 + 5 3 3 0 2 

1.84+ w 

5 21 + - 23 
* 5-Z1 — .18 

338.4 

308 

Fe xvii 

16.780 

0 

3.01 43 7 2 0 

1 m +.59 

A - UA -.oi 

2.45 +\ 3 7 

174.9 

308 



0.5 

4 12 + - 93 
^• az -.90 

1.40+- 22 

2.40+1 

180.7 

308 



1 

c 45+1.46 

1-43*1 

2.47+1 

185.6 

308 



2 

o 77 + 3.39 
°- 7/ -2.60 

1.58- 12 

2.41+1 

193.2 

308 



5 

98 qq+ 15.55 
zo ’^-ll .96 

1-64*1 

2.3S +- 23 

206.3 

308 

0 viii Lya? 

18.969 

0 

3.00 +- 34 

1 . 22+1 

3.70+- 29 

150.9 

130 



0.5 

4 26 +118 

^• zo -l.l 6 

1-43*1 

3.70+- 28 

152.9 

130 



1 

5.22-217, 

1-57*1 

3-70+ ,3 2 7 

155.1 

130 



2 

7-83 43 2 9 9 2 3 

l- 66 *io 

3.62+1 

158.3 

130 



5 

21 - 10 41 , 5 o 9 5 2 o 

| 74+-23 
1 • 7 ^—.20 

3.57 +- 42 

165.0 

130 





Table 5 





Wind Profile Model with Anisotropic Porosity: Fit Results to the XMM-Newton RGS Spectra 



Ion 

X 

hoo 

T* 

Ro 

Normalization 

x 2 

Abins 


(A) 

(«*) 


(R*) 

( 10 -4 photons cm -2 s _1 ) 



Nex Lya? 

12.1339 

0 

1-81*1 

1-61* l 

9 10 +08 
* 5,AU -.08 

24.0 

19 



0.5 

3.92+' 83 

L35+' 2 , 1 , 

o i O+.08 

d.io_ os 

27.3 

19 



1 

5 71+1-23 
7 - 1.00 

1 4T Kl5 
a -^ a -.14 

7 77 +.O 8 
■3- zz — .09 

34.2 

19 



2 

12.61 

1.46 

3.31 

47.2 

19 



5 

64.82 

1.35 

3.44 

66.9 

19 

Fe xvii 

15.014 

0 

1.77 

1.57 

6.39 

129.5 

52 



0.5 

3.29 

1.44 

6.47 

180.4 

52 



1 

4.66 

1.50 

6.52 

249.6 

52 



2 

8.24 

1.57 

6.56 

370.1 

52 



5 

35.96 

1.54 

6.66 

577.2 

52 

Fe xvii 

16.780 

0 

3-38^45 

1-54*39 

9 01 + 07 
:> - ua -.06 

30.3 

32 



0.5 

4.0 1 - 83 

1-99+1 

3 20 + - 07 

^• zu _.07 

42.2 

32 



1 

4-17!. 1 ,! 

2 -24*l 

-2 7Q+-08 

48.0 

32 



2 

4 47+1.87 

^ 7 -1.64 

2-49+1 

3 41 + 09 
^■^ a _.08 

54.5 

32 



5 

4.28 

2.89 

3.59 

62.3 

32 

0 vm Lya 

18.969 

0 

3.14 

1.01 

4.66 

204.3 

72 



0.5 

5.30 

1.35 

4.76 

275.2 

72 



1 

8.14 

1.42 

4.81 

355.1 

72 



2 

16.81 

1.46 

4.87 

478.8 

72 



5 

73.81 

1.37 

4.92 

645.8 

72 

N vii Ly/5 

20.910 

0 

4 Q7+-66 
^■5_1.03 

1-41*1 

1 - 66 +% 

40.5 

36 



0.5 

7 96 +2 - 20 

1-60+1 

1-69+1 

52.9 

36 



1 

12 - 95 -3°34 

1.70+1 

1 75 +11 
1 ' 7 —.12 

62.4 

36 



2 

34.11 

1.62 

1.86 

73.5 

36 



5 

> 100 

1.73 

1.91 

84.4 

36 
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Figure 4. The Chandra MEG (left column) and XMM-Newton RGS (right column) measurements of the Fexvn line at 15.0 14 A; identical to the data shown in 
Figures 2 and 3, but with best-fit isotropic porosity models superimposed. 

(A color version of this figure is available in the online journal.) 
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Figure 5. The 68.3%, 95.4%, and 99.7% confidence limit contours in r*-/zoo parameter space for the anisotropic (left panel) and isotropic (right panel) porosity model 
fits to the Fe xvn line at 15.014 A in the Chandra data. The best fit model is denoted by the filled circle. 




x x 

Figure 6. Comparison of the best fit models for the Fe xvn 15.014 A line in the Chandra data. The model profile flux is shown as a function of the scaled wavelength 
x = (8X/Xo)(c/Voo). Anisotropic porosity is shown in the left panel, while isotropic is shown in the right panel. Note that the best fit anisotropic porosity models have 
a very different shape from the best fit nonporous model because of the enhancement near line center from the “Venetian blind effect.” 


models with large porosity lengths. This is at odds with 
the conclusions of Oskinova et al. (2006), who report that 
anisotropic porosity models provide good fits to the line profiles 
of several bright O stars observed with Chandra , while non- 
porous models do not fit their data as well. 

We show below that the difference in our findings has its 
root in a less obvious assumption of both our models and those 
of Oskinova et al.: the radial upper bound for X-ray emission, 
7? max- Oskinova et al. assume R max = 5* for £ Pup, and as high 
as Rmax = 9/?* for other stars in their sample. On the other hand, 
in this paper we have until now assumed R max = oo. 

5.1. The Effect of the Radial Upper Bound 

In Figure 7 we show a comparison of models illustrating the 
effect of a finite upper radial cutoff to the X-ray emission in 
line profile models. For a given value of r*, a model with a 
finite radial cutoff can produce a much more asymmetric line 
profile. Of course, porosity can mitigate this to some extent. 
However, note that the models with a finite radial cutoff to 
X-ray emission have no flux at all in the wings of the profile, 
since there is no X-ray emitting plasma far out in the wind where 
the highest velocities are. This important difference applies to 
porous models just as much as non-porous models. 


5.2. Comparison with Observed Line Profiles 

To attempt to reconcile our results with those of Oskinova 
et al. (2006), we have performed fits to the Ovm line in the 
RGS spectrum of f Pup assuming the same upper and lower 
radial bounds to X-ray emission that they used: Ro — 1.5 R., and 
R max = 5 R„. We have also fit the same model to the Chandra 
spectra and get comparable results to the RGS fits, but with 
lower statistical significance, and we therefore focus on the 
RGS results in this section. 

We fit a non-porous model as well as models with both 
isotropic and anisotropic porosity. The results of our fits are 
given in Table 8 and the best fit models are shown in Figure 8. 
Because the fits are formally unacceptable, we do not report 
model parameter confidence intervals. Note that the model with 
isotropic clumps preferred no porosity effect in the fit, so we do 
not list it separately. 

The results of our fits assuming finite R max are qualitatively 
similar to Oskinova et al. We find that under these assumptions, 
no value of r* can provide a good fit to the data in the 
absence of porosity, while on the other hand we find a much 
improved fit from an anisotropic porosity model. Nevertheless, 
the anisotropic porosity model with finite upper radial cutoff 
does not provide nearly as good a fit as the nonporous model 
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Table 6 

Wind Profile Model with Isotropic Porosity: Fit Results to the Chandra HETGS Spectra 


Ion 

(A) 

hoo 

( R *) 

T* 

Ro 

(*,) 

Normalization 
(10~ 4 photons cm -2 s -1 ) 

C 

Abins 

Mg xn Ly a 

8.4210 

0 

1 22 +-78 
I.ZZ_ 44 

1-342% 

0.2942° 0 % 

186.5 

188 



0.5 

1-312% 

1 -382.18 

0.3002°o 1 2 9 9 

186.8 

188 



1 

1.462 58 

i .392.15 

0.3002°o 1 2 9 8 

187.1 

188 



2 

1 70+1-42 
1 - /y -.80 

1.412% 

0.3002°o 2 2 ° 8 

187.7 

188 



5 

T i q+4.33 
1 — 1.84 

1 -472.10 

0.3012°o 2 2 ° 8 

189.7 

188 

Nex Ly a 

12.1339 

0 

2 . 012 % 

1 4C+.13 
A - q °— .08 

2-712.09 

191.4 

176 



0.5 

2.392% 

1 40 +IO 
A — .08 

2.692% 

191.8 

176 



1 

2-802% 

1 51 +09 
AJA - .07 

2 69 +10 
z - oy -.08 

192.5 

176 



2 

3.862% 

1 54+-07 
a - d ^-.08 

2 6 Q +1 ° 
z - ov -.08 

195.0 

176 



5 

q 94+2.82 

1-612% 

2.692% 

208.5 

176 

Fe xvii 

15.014 

0 

1.942 ' 33 

1-552.12 

5.242% 

280.8 

308 



0.5 

2 . 182 4 3 2 9 

1.582.io 

5.242% 

282.0 

308 



1 

2.462% 

1.612.13 

5.242% 

283.4 

308 



2 

3-092.66 

1 63 +08 

1 .OO-09 

5-242% 

286.4 

308 



5 

5.702 2 1 % 

1 74+.07 
A,/ ^-.09 

5.232% 

297.2 

308 

Fe xvii 

16.780 

0 

3.012% 

1 01 + - 59 
A - UA -.oi 

2.452% 

174.9 

308 



0.5 

3.592% 

L352.3 8 ! 

2.392% 

174.1 

308 



1 

4-3B2V 1 o 1 6 

1 .422% 

2.432.15 

173.3 

308 



2 

6 52 +1 " 

o.oz—i ? 5 

1 .482.il 

2.402% 

172.6 

308 



5 

17-882 6 6 5 oo 

L602.ii 

2.432.14 

176.3 

308 

O viii Lya 

18.969 

0 

3.00+% 

1.22+% 

3-702.35 

150.9 

130 



0.5 

7J oc+1.05 
3.5D_1 10 

1 22 + - 38 
A " r ’ z -.20 

3.702% 

150.9 

130 



1 

4 42 +1 - 74 
— 1.47 

1 4S+- 34 

3-712.29 

151.0 

130 



2 

c q^+ 2.88 
J . oD _2 ig 

1.562% 

3.622 % 

151.1 

130 



5 

1 "3 22+8-74 
.80 

1.692% 

3.602% 

152.1 

130 



X 


Figure 7. Comparison of three model profiles: a model with no porosity and no 
upper radial cutoff; a model with no porosity and upper radial cutoff of 5 /?*; 
and a model with anisotropic porosity, h 0 0 = R *, and an upper radial cutoff of 
5 R *. All of the models assume r* = 3. Note that both of the models with a 
cutoff have corresponding cutoffs to X-ray emission in the wings of the profile. 

of Table 3 with no upper radial cutoff, which is also shown in 
the first panel of Figure 8 for comparison. The poor fit of both 
models with 7? max = 5 7?* to the data is due in large part to the 
lack of flux in the wings of the profile, especially the red wing. 
This lack of flux in the wings is a necessary consequence of a 
finite radial upper cutoff. 


5.3. Constraints on the Radial Upper Bound 
from Forbidden Line Strength 

An obvious way to obtain a strong constraint on the presence 
of hot plasma at large radii is by modeling the weak forbidden 
lines of the He-like triplets, which are sensitive to the radial 
distribution of X-ray emitting plasma (Kahn et al. 2001). This 
sensitivity comes about because the upper level of the forbidden 
line is metastable, and can be depopulated by photoexcitation 
from the star’s intense UV radiation field, the strength of which 
is radially dependent due to geometrical dilution. 

The forbidden-to-intercombination line ratio is predicted for 
different He-like ions as a function of radius in Figure 3 of 
Leutenegger et al. (2006). It is evident from this plot that the 
forbidden lines of N vi and O vn are strongly suppressed for 
r < 207?*, which already suggests that an upper radial cutoff 
of 7? max = 57?* is likely not compatible with the detection of 
weak forbidden lines. However, as shown in Leutenegger et al. 
(2006), it is straightforward to include the radial dependence of 
the // i ratio in the integrand of the line profile calculation in 
Equation (1) and test the effects of an upper radial cutoff against 
the measured He-like complex data. 

We fit this model to the O vii and N vi triplets in the RGS 
spectrum of £ Pup. Resonance scattering was included in our 
models, as described in Leutenegger et al. (2007). The weak 
Cvi Ly/3 line (which falls on the red wing of the resonance 
line of N vi) was included in our fit to the N vi complex, as in 
Leutenegger et al. (2007). Both complexes were fit under the 
competing assumptions that 7? max = oo or 207?*. In order to 
make a more direct comparison with the radial distribution of 
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Table 7 

Wind Profile Model with Isotropic Porosity: Fit Results to the XMM-Newton RGS Spectra 


Ion 

J, 

(A) 

hoo 

( R *) 

T* 

Ro 

(R*) 

Normalization 
( 10 -4 photons cm -2 s _1 ) 

X 2 

A^bins 

Nex Lya? 

12.1339 

0 

1 . 8 Pf 2 

l-61-\ 5 8 

0 iq +.08 
■ 5,au -.08 

24.0 

19 



0.5 

2.07jj 27 

1 61- ‘,6 

3 10 +08 
- 3 - 1 -.08 

24.0 

19 



1 

2 41+- 45 

i-«--i s 5 

3 10 +09 
■ 5,au -.08 

24.1 

19 



2 

3-42^i 

1 - 60 - 1 l 4 3 

3 1 1 + 08 
- 5 - 1 a -.09 

24.2 

19 



5 

1 D QQ+6.62 

1.58-12 

3 16 +08 
■ 3 - ad -.08 

26.3 

19 

Fe xvii 

15.014 

0 

1.77 

1.57 

6.39 

129.5 

52 



0.5 

2.02 

1.58 

6.40 

131.4 

52 



1 

2.33 

1.59 

6.40 

133.6 

52 



2 

3.23 

1.60 

6.40 

139.3 

52 



5 

9.96 

1.61 

6.44 

175.6 

52 

Fe xvii 

16.780 

0 

3.38^ 

1.54-33 

3 01 + 07 
- 5 - ua -.06 

30.3 

32 



0.5 

3.94^1 

1.65-1 

3 01 + 07 
- 5 - ua -.06 

31.1 

32 



1 

4.50— 8 g j 

1 73 + - 23 

3.03-°l 

32.2 

32 



2 

5-61^7 

1 . 86-1 

3.06-° 0 7 7 

34.4 

32 



5 

8.61^ 3 7 0 6 

2-13-.1 

3 19 +07 

40.4 

32 

O viii Lya? 

18.969 

0 

3.14 

1.01 

4.66 

204.3 

72 



0.5 

3.64 

1.57 

4.64 

198.0 

72 



1 

4.35 

1.51 

4.64 

193.4 

72 



2 

6.38 

1.57 

4.63 

188.9 

72 



5 

20.20 

1.65 

4.67 

222.3 

72 

N vii Ly/1 

20.910 

0 

4.93 -“ 3 

1 41 +.62 

a -^ a -.40 

1 . 66 -% 

40.5 

36 



0.5 

6-18-i.4 9 6 

1.65- 4 5 6 i 

1 . 66-1 

40.6 

36 



1 

7 - 42 -l 6 9 7 l 

1 7Q+.36 
1 ‘ ''—.36 

1 . 66-1 

41.0 

36 



2 

10.16L 4 2 2 7 6 6 

1 Q5+-33 
A - y:) -.30 

1 . 66-1 

42.4 

36 



5 

26-09+1 

2 14 + - 38 
z ' a ^-.29 

1.69-1 

49.0 

36 


Table 8 

Effect of Radial Bounds of Rq = 1 .5 R* and R max = 5R* on Fit to O viii Lya? 
in XMM-Newton RGS Spectrum 


hoo 

(■ R *) 

r* 

Type 

Normalization 
(10~ 4 photons cm -2 s -1 ) 

x 2 

M)ins 

0 

2.34 

np, iso 

4.404 

1421.0 

72 

1.03 

18.46 

aniso 

4.506 

769.6 

72 


X-ray emitting plasma assumed in Oskinova et al., we also tried 
to fit the data with models using their values of Rq = 1.5 R* and 
J? max = 5R :t . and including the effects of porosity by allowing 
hoo to be a free parameter. The results of our fits are given in 
Table 9, and the models are shown in Figures 9 and 10. 

The models with a finite radial cutoff clearly do not 
produce enough emission from the forbidden line to fit the data. 
On the other hand, the models with no radial cutoff fits the // i 
line ratios (and the line complexes) quite well. We thus conclude 


that a R max = 5 R* cutoff is strongly disfavored by the observed 
forbidden line strengths as well as the wings of individual lines. 

6. DISCUSSION 

One of the most important results of our emission line profile 
modeling is that anisotropic porosity models are disfavored. 
This can be understood in terms of the markedly different 
profile shapes predicted by anisotropic porosity models, with 
a strong bump at line center due to the “Venetian blind effect,” 
as exemplified in Figure 6. There is no such bump evident in 
any of the profiles we have modeled, and this is reflected in our 
formal statistical constraints. 

On the other hand, isotropic porosity models can provide 
good fits to the lines we studied, provided the porosity lengths 
are no larger than R f . Porosity length and mass-loss rate are 
degenerate to a certain extent, since they both influence the 
degree of profile asymmetry. We have quantified this degeneracy 
for one line in Figure 5. Large porosity lengths (h 0 0 > R. f .) and 
their associated large adjustments to derived mass-loss rates 


Table 9 

He-like Triplet Models with Finite Radial Upper Cutoff: Fit Results to the XMM-Newton RGS Spectra 


Ion 

r* 

Ro 

(R*) 

^max 

(■ R *) 

hoo 

(R*) 

7 0 ,* 

G 

Normalization 
(10 -4 photons cm -2 s -1 ) 

x 2 

^bins 

Nvi 

4.79 

2.28 

00 

0 

00 

1.162 

17.52 

181.5 

142 


5.70 

4.0 

20 

0 

00 

1.176 

16.96 

381.5 

142 


26.28 

1.5 

5 

0.72 

00 

1.085 

16.05 

1435.9 

142 

O vii 

3.96 

1.59 

00 

0 

2.68 

1.005 

7.582 

146.8 

109 


7.75 

4.0 

20 

0 

8.32 

1.174 

7.408 

311.1 

109 


25.57 

1.5 

5 

0.77 

0.85 

1.095 

7.171 

740.8 

109 


12 


The Astrophysical Journal, 770:80 (17pp), 2013 June 10 


Leutenegger et al. 



18.8 18.9 19.0 19.1 

Wavelength (A) 




Figure 8. Three model fits to the Oviii Lya line in the XMM-Newton 
RGS spectrum of £ Pup. RGS2 data are shown, but the fits are performed 
simultaneously to data from both RGS. The first panel shows the nonporous 
fit reported in Section 4. The second panel shows a non-porous model with 
Rq = 1.5 R* and R max = 5 R*. The third panel shows a model with anisotropic 
clumps, and with the same radial cutoffs. The last two panels correspond 
approximately to the models shown in Figure 8 of Oskinova et al. (2006). 

(A color version of this figure is available in the online journal.) 

are disfavored by the data. Ruling out large porosity lengths 
is possible because porosity length and mass-loss rate are not 
fully degenerate parameters, leading to measurable differences 
in profile shape, as shown in Figure 6. Smaller non-zero porosity 
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Wavelength (A) 
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Wavelength (A) 
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Figure 9. XMM-Newton RGS spectrum of the O vn complex of £ Pup, together 
with best fit models assuming R max = oo (first panel), R md \ = 20 R* (second 
panel), and R max = 5 R* (third panel). The three vertical dashed lines indicate 
the positions of the rest wavelengths of the resonance, intercombination, and 
forbidden lines from left to right. The model shown in third panel includes 
the effects of anisotropic porosity, and is comparable to the model favored 
in Oskinova et al. (2006). Both models with finite radial upper cutoffs 
underpredict emission from the forbidden line, which is primarily formed at large 
radii. 

(A color version of this figure is available in the online journal.) 

lengths are allowed, but lead to only modest adjustments to 
derived mass-loss rates. 

We have compared our results to those of Oskinova et al. 
(2006), who also compare both porous and non-porous models 
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Figure 10. XMM-Newton RGS spectrum of the N vi complex of £ Pup, together 
with best fit models assuming /? max = oo (first panel), R max = 20 R* (second 
panel), and R max = 5R* (third panel). The three vertical dashed lines indicate 
the positions of the rest wavelengths of the resonance, intercombination, and 
forbidden lines, from left to right. The model shown in third panel includes 
the effects of anisotropic porosity, and is comparable to the model favored in 
Oskinova et al. (2006). Both models with finite radial upper cutoffs underpredict 
emission from the forbidden line, which is primarily formed at large radii. The 
weak feature at 28.4656 A is C vi Ly/J. 

(A color version of this figure is available in the online journal.) 


to X-ray spectra of £ Pup, but reach the opposite conclusion that 
models assuming anisotropic porosity are allowed, and non- 
porous models disfavored. We showed that this difference in 
conclusions is a direct result of a difference in assumptions 
regarding the spatial distribution of X-ray emitting plasma. 
Specifically, we assumed in our modeling that X-ray emission 
extends to effectively infinite radius, while Oskinova et al. 
assumed that X-ray emission turns off above 5 /?* . 

In Figure 7 we compared the effect of different assumptions 
regarding /? max , and found that models with a finite radial cutoff 
have much less flux in the wings of the line profile than models 
with no cutoff, which is a straightforward consequence of the 
fact that all emission at large redshift and blueshift originates 
at large radii. We fit models with different assumed values for 
Rmax to the O viii Lya line, and found that models with a finite 
radial cutoff do not fit the wings of the observed line profile of 
? Pup. 

We further probed the radial distribution of X-ray emission 
in the wind of £ Pup by fitting the He-like triplet lines of O vii 
and N vi with models accounting for the radial dependence of 
the forbidden-to-intercombination line ratio (Leutenegger et al. 
2006). As is evident from Figures 9 and 10, models with a 
finite radial cutoff do not reproduce the observed forbidden line 
flux, which indicates that there is significant X-ray emission 
originating from R > 20/?*. 

The finding of significant X-ray emission from large radii 
is one of the major, if unintended, results of this study. One 
might be surprised by this result based on the simple argument 
that X-rays should originate where the wind is accelerating 
and can generate strong shocks. Runacres & Owocki (2002) 
have performed hydrodynamic simulations specifically aimed 
at understanding the outer wind structure of O stars, and they 
find that X-ray emitting shocks can be generated at tens of stellar 
radii. Furthermore, the cooling time for few MK plasma formed 
at r > 10/?* in the wind of t, Pup is comparable to the flow 
time, so even if new shocks are not generated at very large radii. 
X-ray emission can still persist (Feldmeier et al. 1997a). 

Under similar modeling assumptions as Oskinova et al. 
(2006), in particular that X-ray emission is cut off above 5 /?*, we 
too find that porous models with anisotropic clumps are favored 
over non-porous models. But we also show that such a small 
cutoff radius for X-ray emission is not consistent with either 
the observed flux in the wings of X-ray emission lines, or the 
strength of forbidden line in He-like triplets. These constraints 
are well matched by models without an arbitrary X-ray cutoff. 
With such extended emission, fitting line profiles favors models 
with either no or a modest isotropic porosity, with anisotropic 
porosity strongly disfavored. 

Cohen et al. (2010) have derived a mass-loss rate for £ Pup 
from modeling of X-ray line profiles under the assumption that 
porosity effects are negligible. Taking the conclusions of the 
present work at face value, the correction to this mass-loss rate 
is at most 40% if moderate porosity effects are present (/ioo ~ 
/?*). This conclusion is further supported by the fact that the 
uncorrected X-ray mass-loss rate of Cohen et al. (2010) agrees 
with mass-loss rates from non-X-ray observational diagnostics 
within uncertainties (Najarro et al. 2011; Bouret et al. 2012). 

We thus conclude that, at least for f Pup, X-ray line profiles 
are a good independent diagnostic of mass-loss rates, and that 
they are not subject to strong systematic errors from clumping 
on any scale. The other O stars observed by Chandra and XMM- 
Newton should ideally also be subjected to a similar study, 
although it would perhaps be surprising to find major differences 
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in porosity effects in the wind of £ Pup in comparison with 
other O stars. In fact, early O supergiants like £ Pup with 
their high mass-loss-rate winds are the stars for which porosity 
effects are expected to be the strongest. Because the essence of 
porosity is the optical thickness of individual clumps, O stars 
with lower mass-loss rates than £ Pup require more extreme 
clump properties just to produce the same porosity effect. 

The relative robustness against wind inhomogeneity effects 
of X-ray line profile measurements of mass-loss rates suggests 
their use as a primary mass-loss rate diagnostic, on an equal 
footing with traditional diagnostics such as Ha, thermal ra- 
dio emission, and UV absorption line profiles. However, X-ray 
line profile measurements are currently only possible for the 
brightest stars in the Galaxy. Until the advent of much larger 
X-ray spectroscopic observatories, with square meter effective 
areas and which could undertake a large-scale survey of O star 
X-ray spectra, we suggest a program of benchmarking as many 
stars as possible with all available mass-loss rate diagnostics. 
Only by continuing the tradition of a multiwavelength, multidi- 
agnostic approach can we hope to disentangle the observational 
signatures of mass-loss and wind structure. 
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APPENDIX 

EXPANDED FIT RESULTS 

In this section in Figures 1 1-20 we show plots of all the fits 
reported in Tables 4-7. For each line that we fit, we show plots of 
anisotropic porosity models in one figure, and isotropic porosity 
models in another figure. We repeat the non-porous model in 
both the anisotropic and isotropic porosity fit figures. To keep 
the figures legible, we only show models with h 0 0 = 0, 1 , and 
5 R*. 



Figure 11. The Chandra MEG measurements of the Mg xn line at 8.421 A, with 
the best-fit anisotropic porosity models superimposed. The red, green, and blue 
models assume h 0 0 =0, 1, and 5 R *, respectively. 

(A color version of this figure is available in the online journal.) 



Figure 12. The Chandra MEG measurements of the Mgxn line at 8.421 A, 
with the best-fit isotropic porosity models superimposed. The red, green, and 
blue models assume h 0 0 = 0, 1, and 5 R *, respectively. 

(A color version of this figure is available in the online journal.) 
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Figure 13. The Chandra MEG (upper panel) and XMM-Newton RGS (lower 
panel) measurements of the Nex line at 12.134 A, with the best-fit anisotropic 
porosity models superimposed. The red, green, and blue models assume h OQ = 
0, 1, and 5 R *, respectively. 

(A color version of this figure is available in the online journal.) 
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Figure 14. The Chandra MEG (upper panel) and XMM-Newton RGS (lower 
panel) measurements of the Nex line at 12. 134 A, with the best-fit isotropic 
porosity models superimposed. The red, green, and blue models assume h OQ = 
0, 1 , and 5 R * , respectively. 

(A color version of this figure is available in the online journal.) 
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Figure 16. The Chandra MEG measurements of the Fexvn line at 16.780 A, 
with the best-fit isotropic porosity models superimposed. The red, green, and 
blue models assume h 0 0 = 0, 1, and 5 R *, respectively. 

(A color version of this figure is available in the online journal.) 
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Figure 15. The Chandra MEG measurements of the Fexvn line at 16.780 A, 
with the best-fit anisotropic porosity models superimposed. The red, green, and 
blue models assume h^ = 0, 1, and 5 R *, respectively. 

(A color version of this figure is available in the online journal.) 



18.8 18.9 19.0 19.1 

Wavelength (A) 


Figure 17. The Chandra MEG (upper panel) and XMM-Newton RGS (lower 
panel) measurements of the O vm line at 18.969 A, with the best-fit anisotropic 
porosity models superimposed. The red, green, and blue models assume h Q Q = 
0, 1, and 5 /?*, respectively. 

(A color version of this figure is available in the online journal.) 
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Figure 18. The Chandra MEG (upper panel) and XMM-Newton RGS (lower 
panel) measurements of the Ovm line at 18.969 A, with the best-fit isotropic 
porosity models superimposed. The red, green, and blue models assume hoc = 
0, 1, and 5 /?*, respectively. 

(A color version of this figure is available in the online journal.) 



Figure 19. The XMM-Newton RGS measurements of the N vn line at 20.9 1 0 A, 
with the best-fit anisotropic porosity models superimposed. The red, green, and 
blue models assume hao = 0, 1, and 5 /(*. respectively. 

(A color version of this figure is available in the online journal.) 



Figure 20. The XMM-Newton RGS measurements of the N vn line at 20.910 A, 
with the best-fit isotropic porosity models superimposed. The red, green, and 
blue models assume /too = 0, 1, and 5 R t , respectively. 

(A color version of this figure is available in the online journal.) 
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